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Abstract 

This paper studies cooperative spectrum sensing in cognitive radio networks where secondary users 
collect local energy statistics and report their findings to a secondary base station, i.e., a fusion center. 
First, the average error probability is quantitively analyzed to capture the dynamic nature of both 
observation and fusion channels, assuming fixed amplifier gains for relaying local statistics to the 
fusion center. Second, the system level overhead of cooperative spectrum sensing is addressed by 
considering both the local processing cost and the transmission cost. Local processing cost incorporates 
the overhead of sample collection and energy calculation that must be conducted by each secondary 
user; the transmission cost accounts for the overhead of forwarding the energy statistic computed at each 
secondary user to the fusion center. Results show that when jointly designing the number of collected 
energy samples and transmission amplifier gains, only one secondary user needs to be actively engaged 
in spectrum sensing. Furthermore, when number of energy samples or amplifier gains are fixed, closed 
form expressions for optimal solutions are derived and a generalized water-filling algorithm is provided. 



I. Introduction 

To alleviate inefficient allocation of radio frequency (RF) spectrum, cognitive radios have recently 
been proposed to coexist with primary (or licensed) users of spectral bands while not causing harmful 
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interference HI ill- Current proposals for secondary networks require cognitive users to conduct spectrum 
sensing so that they can detect unused spectral bands and avoid interfering with a primary system. To 
improve detection reliability in fading conditions, multiple secondary users can cooperate in spectrum 
sensing and take advantage of spatial diversity |[3lHl . 

In secondary networks where users communicate with a local secondary base station as illustrated in 
Fig. [T] the system level performance and design of cooperative spectrum sensing must 1) account for the 
dynamic nature of both the observation and fusion channels, i.e., the channel between the secondary and 
primary users and the channel between the secondary user and the secondary base station, respectively; 
and 2) balance the gains offered by spectrum sensing against its computational and transmission costs. 
In this paper, we address both these concerns in evaluating and designing spectrum sensing schemes for 
secondary networks. 
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Fig. 1. Topology of cooperative spectrum sensing in cognitive radio networks. 

In (H, a logic OR fusion rule for hard-decision combining was presented to cooperatively detect the 
primary user. The AF cooperative strategy was used in J6[ to improve spectrum agility and allow two 
secondary users to communicate with each other. An optimal linear detector for cooperative spectrum 
sensing was proposed in [3), where the received signals at the fusion center were optimally weighted 
for global fusion. In Q, a linear quadratic fusion rule based on a detection criterion was proposed for 
spectrum sensing by modeling received signals as correlated log-normal random variables. Based on 
our knowledge, these and other prior studies do not focus on system-level performance of cooperative 
spectrum sensing that accounts for the dynamic nature of both the observation and fusion channels. 

Low-energy overhead cooperative spectrum sensing was studied in (8). Optimally allocated powers were 
computed without taking into account the underlying system level cost of sensing. Our work on energy- 
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constrained spectrum sensing is motivated by (9), where detection problems accounted for constraints 
on expected cost due to transmission and measurement. We build on these formulations here to design 
energy-constrained cooperative spectrum sensing. 

In our system model, secondary users forward local energy statistics to a secondary base station 
using amplify-and-forward (AF) over parallel access channels. We first address the impact of dynamic 
observation and fusion channels by analyzing the average error probability for cooperative spectrum 
sensing considering both additive white Gaussian noise (AWGN) and Rayleigh fading conditions. Results 
show that detection performance can be maintained in the low and moderate fusion signal-to-noise ratio 
(SNR) regimes when fusion channels are reliable, whereas fading on the secondary users' observation 
channels provide spatial diversity. 

Next, we address the system level energy cost of sensing by considering two major factors: Local 
processing cost due to sample collection and local energy calculation and transmission cost due to 
forwarding local statistics to the fusion center. We present two optimization problems to find the number 
of energy samples that must be collected at each secondary user and the appropriate amplifier gain that 
each secondary user must use for AF relaying of the local energy statistic. When jointly optimizing both 
the number of samples and amplifier gains, we show that only one secondary user must be actively 
engaged in spectrum sensing. When either the amplifier gains or the number of samples is fixed, we find 
closed-form optimal solutions and propose a generalized water-filling approach to energy-constrained 
cooperative spectrum sensing. 

The remainder of the paper is organized as follows: Section II describes our system model. Section 

III presents the average error probability for various observation and fusion channel conditions. Sections 

IV and V collectively present our results for energy-constrained spectrum sensing: Section IV addresses 
the optimization for minimization of global error probability while Section V provides the optimization 
for minimization of system level cost. Simulation results are presented in Section VI and we conclude 
the paper in Section VII. 

In this paper, we use the following notation: column vectors are denoted by boldface lowercase letters, 
i.e., x = [x\,X2, • • • , x n ] T and Xi is the ith entry of x. = [0, 0, • • • , 0] T and 1 = [1, 1, • • • , 1] T . I is the 
identity matrix. (-) T and (•)' denote the transpose and conjugate transpose operation, respectively. 
denotes the £2 norm of x. x y denotes the generalized inequality, i.e., Xi > 0. Z 1 } and W) denote the 
set of nonnegative integer and real n-vectors, respectively. |<S| denotes the cardinality of a set S. \-~\ and 
[•J denote the ceiling and floor operations, respectively. 
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II. System Model 

A. Communication Model 

We consider a network model in Fig. [T] where secondary user conducts local spectrum sensing and 
transmits its local energy statistic to the fusion center using AF on parallel access channels (PAC). The 
received signal for secondary user i at the fusion center is shown in Fig. |2j i.e., 



Vi = QihiXi + Vi 



(1) 



where x^ is the energy of received signal at the secondary user i; gi is the amplifier gain for the secondary 
user i; hi is the channel gain between secondary user i and the fusion center and v% is independent and 
identically distributed (i.i.d.) Gaussian noise, i.e., v\ ~ CJ\f(0,a%) and is independent of Xj. We assume 
that hi is known at the fusion center (e.g., via channel estimation) and remains constant during the sensing 
period. We can then rewrite ([]]) in a matrix form as 



y = Hx + v, 



where H = diag{gxhi, g 2 h 2 , ■ ■ ■ ,g n h n }- 



(2) 
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Fig. 2. Cooperative spectrum sensing in cognitive radio networks. 



(3) 



B. Local Energy Statistic 

For secondary user i, (1 < i < n), the hypothesis test for Xj is given as 

H : x i = (l/K i )Y,kUWi(k)\ 2 
%i : Xi = (l/«») YJk=i \his(k) + rn(k)\ 2 , 
where m is the number of samples, s(k) is the transmitted signal from the primary user and rii(k) is the 
noise received by secondary user i. We assume s(k) is complex PSK modulated and i.i.d. with mean zero 
and variance a 2 ; hi is the channel gain between the primary user and secondary user i and is assumed to 
be constant during the cooperative spectrum sensing period; and n,i(k) is i.i.d. Gaussian noise with mean 
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zero and variance a\ and is independent of s(k). We define the local received SNR at the secondary user 
i as ji = a1\hi\ 2 / a^. When m is large, Xi can be approximated as Gaussian random variable Q, i.e., 

We assume here the local received SNR 7,- is known at secondary user z. In IEEE 802.22, 7, can be 
estimated from pilot signals periodically transmitted by primary users ifTOl 1 , 



Given this system model, we see that £3 = E{x?} = [1 + l/re, + tt\ (7$ + 2 (1 + 7i] a„, where 

7ro = P(7^o) an( l TTi = P(^i) are the probabilities that spectrum is idle and occupied, respectively. In 
cognitive radio networks, the received primary user power measured by the secondary user is expected 
to be very small ifTTI . i.e., ji <C 1. Additionally, the number of samples is expected to be more than 
a few, i.e., Ki S> 1. Thus, we can approximate the transmitted power for the secondary user i as V% = 

C. Optimal Fusion Rule 

Under hypothesis T~Lq and Hi, the received signal y has a Gaussian distribution, i.e., 

Ho : y ~ N {Hlal S ) 
Hi : y ~JV(H(1 + 7)o* S x ) , 

where S = HSHt^ + cr^I and Si = HS(I + 2r)Ht^ + a*l, here, T = diag{ 7 i, 72, • • • , 7„} and 
S = diag{l/«i, 1/^2, • • • , 1/K n }. Without loss of generality, assume that -kq = m = 0.5. Then, optimal 
(maximum a posteriori probability) likelihood ratio test (LRT) is given as: 

log « 1'0. © 
p(y|^o) % 

Since 7$ <C 1 and 3> 1, then, 7f/ACj ~ and we have S ~ Si. Thus, the optimal LRT can be 
approximated as 

T{y) = (H 7 )tS : 1 y t r, (7) 

Ho 

where r = (H7)tS " 1 H(l + 0.57)<t^. Furthermore, we note that E{T(y)|% } = (H/yJ+Ep ^lo*. 
E{T(y)|^i} = (H 7 )tS- 1 H(l + 7 ) C 72 and Var{T(y)|H } = Var{T(</)|^i} = (H 7 )tS- 1 H 7 . 



'Uncertainty in the knowledge of local received SNR would affect the design of cooperative spectrum sensing. We will 
investigate this important issue in the future. 
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With this preparation, it can be shown that the error probability is given aj^| 

where a% = o~^/a^ and Q(x) = -7= J x °° exp(— t 2 /2)di. It is also easy to see that the asymptotic error 
probability expressions when the number of samples or amplifier gains approach infinity are given by 



and 




respectively 

D. System Level Cost for Cooperative Spectrum Sensing 

In this paper, we consider system level cost for cooperative spectrum sensing in cognitive radio 
networks. This system level cost has contributions from three components: Local processing; transmission; 
and reporting and broadcasting. 

• Local processing cost includes the energy consumed by the secondary user in receiver RF scanning 
and local energy calculation. For simplicity, we assume that the local processing cost C p i(-) for 
secondary user i is a linear function of the number of samples |12|. i.e., C p i{ni) = cqki, where Co 
is the local processing cost per sample. 

• Transmission cost is the transmit power required from a secondary user to transmit the local 
calculated energy to the fusion center. Here, we assume that this cost for secondary user i is given 
as C ti (gi)=Vi = ^gf. 

• For optimal system design, the fusion center needs to know the local received SNR for each secondary 
user. In practice, this means that secondary users will report their local received SNRs to the fusion 
center. The fusion center then determines optimal allocations (number of samples and/or amplifier 
gains) to each secondary user and then broadcasts them to all secondary users. In this paper, we 
assume that this total reporting and broadcasting cost C r b is fixed; thus we do not consider it in the 
optimization problem. 



2 It is worth mentioning that we can reach same optimization formulation by using Neyman-Pearson criterion to maximize 
global detection probability. Here we present the global error probability for the sake of simplicity. 
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The system level cost during the cooperative spectrum sensing (aside from C r b) is given as 

n n n 

C(K,g) = ^Cpi(Ki) + ^C t i{gi) = ^2 {c Ki + Ci9i) ■ 

i=l i=l i=l 

III. Average Error Probability 

We assume in this section that the amplifier gains and the number of samples collected at each secondary 
user are fixed and not adjusted according to the channel gains. We will discuss adapting amplifier gains 
and number of samples in the subsequent sections. From ([8]), we see that the average error probability 
can be calculated as 

Pe,avg = E 7jft {P e | 7ih } . (11) 

To simplify the calculation of the average error probability, we consider the following alternate 
expression for the Q function HI Q(x) = \ £ exp I — 9 x 5 t )d<f), x > 0. When the local received 
SNRs 7j and fusion channel gains | hi \ 2 are independent, respectively, we can simplify the average error 
probability in ([TTI ) as 

i p-k/2 n 

Pe.avg = - / T#i(0)cl0, (12) 

where Bi(4>) = / °°exp (^qr) p yi (s)P\h^(t)dsdt and Ai(s,t) = -| • g |f ■ Here, p 7i (s) and 
p\ h .\2(t) are PDFs of 7^ and |/ij| 2 , respectively. If we further assume that = g, Ki = k, ji and |/ij| 2 
are i.i.d., respectively, i.e., p yi (s) = p 7 (s) and p\ h .\2(t) = p\ h ^(t), then we have £>j(</3) = £>(</>), VI In 
this case, the average error probability in (TTTb reduces to 

Pe,avg = - r /2 [fi(0)] n d0. (13) 
7T Jo 

Based on this, we see that P ej avg is a decreasing function of n, which indicates that in a power uncon- 
strained cognitive radio network, global error performance can be improved by increasing the number 
of secondary users. This statement follows since A(s,t) < and £>(<7j) < / °° / °° p 1 (s)p\h\^(t)dsdt = 1. 
In general, a closed-form expression of P e ,avg is difficult to obtain. However, only elementary functions, 
such as exponential and Q(-), are involved in the integral calculation; the average error probability can 
thus readily be found numerically. 

Remark: To gain more insight, we investigate an upper bound for average error probability. Since 
Q(x) < iexp(-x 2 /2), the upper bound can be obtained as 

1 - 

P e,avg = 2 11-^- < 14 ) 

i=l 



8 



where Mi = J* °° J* °° exp[Ai(s,t)]p li (s)p\h i \2(t)dsdt. Assume for simplicity gi = g, K{ = k, 74 and 
\hi\ 2 are i.i.d., respectively. Then, we have Mi = M, Vi, and 

Pe.avg = ^M". (15) 

It is readily evident that when g — > 0, P e ,avg — ► o- This * s not surprising since when the amplifier gains 
are low, the fusion center will not be able to make a global decision due to the lack of local energy 
statistic. Next, we use ([T3l and (fl3T ) to evaluate the average error probability for cooperative spectrum 
sensing for the three channel scenarios shown in Table HI 

TABLE I 

Three Channel Environments for Performance Evaluation 





Observation channels 


Fusion channels 


Channel Environment I 


AWGN 


Rayleigh fading 


Channel Environment II 


Rayleigh fading 


AWGN 


Channel Environment III 


Rayleigh fading 


Rayleigh fading 



1) Channel Environment I: In this scenario, 7i = 7 and pi/^^i) = exp(— t) since the observation 
channel is AWGN and the fusion channel is exponential Rayleigh fading. After some manipulations, we 
have 



B(4>) = exp 



K7 



) 5111 



K7 



> sin 



r 



where ^1 (a, b) = exp ( — x + -£^)dx, (a, b > 0). After calculating B((f>), we substitute it in (TT31 to 
obtain the average error probability. It is interesting to note that a similar definition of (<p, a, b) can 
be found in (6J. Furthermore, the upper bound is given as 



e,avg 



2 exp 



nivy 



_o ~2 1 n 



When g — > 00, we see that ^Lgidoo) = \ exp 



nn-y- 



*1 



. This indicates that when the fusion channel is 



perfect, average error performance is limited by local observed energy statistic. 

2) Channel Environment II: In this scenario, p 7i (s) = i exp(— 4) and hi = 1. After some manipula- 
tions, (using eq.(3. 322.2) in |[l4l ). we obtain 

B((f)) = V 8nc sin (j> exp (2c sin 2 (p) Q (2y/csm(j>) , 

where c=-^(^ + % J . Furthermore, the upper bound is 



i( 2 ) 



P^avg = \ (Svrc) " /2 exp (2nc) [Q {2yfc ) 



When g -> oo, we see that c -»• l/(«7^) and P^^oo) = 5 I # ) exp ( ^ ) Q ( ^ 



2 



7v^ 



. Again, 



we see that the average error performance is limited by local observed energy statistic when g — > oo. 

3) Channel Environment III: In this scenario, p 7l (s) = iexp(— 4) and p^.^i) = exp(— t). After 
some manipulations, we have 



„. ,, / — /2sin 2 <A /sin 2 ^ cr 2 sin 2 <A 
8(0) = ^exp ) * 2 ) 



where ^ 2 (a, 6) = / °° {a + §) 7 exp (-x + f ) Q f 2 (a + § ) 7 J dx, (a, b > 0). Furthermore, the 
upper bound is 

f>( 3 ) 1 to \n/2 ( 2n \ [ / 

P e ,av g =2(^) -P^J [* 2 (- 



1 6T. 



,K7 2 ' <7 2 7 2 

~(3) ~(2) 

When g — )• oo, we see that P e!avg (<7oo) = Pe.avglfi'oo)- This is primarily due to the fact that when g — > oo, 
the fusion channel no longer impacts the average error performance. 

IV. Optimization: Minimization of Error Probability 

In this section, we aim to minimize the error probability for the system model in Fig. [2] subject to a 
system level cost constraint of sensing. Specifically, we determine the appropriate number of samples and 
amplifier gains for each secondary user and consider the following two scenarios for this optimization 
problem: 

1) Scenario A: First, we consider the system level cost constraint. Hence, the optimization problem 
is formulated as: 

min V e (n,g) 
K,g 

s.t. C(n,g) < C, k G Z™, g G Tl n + , (16) 

where C is the system level cost constraint. Here we denote the optimal solution of (fT6l) as 
{ K *pT > 9pT ) anc ^ ^ e minimum error probability as p^, opt,1 \ 

2) Scenario B: In some applications, local sample collection for each secondary user may be scheduled 
in a fixed time slot. This indicates the number of samples is upper bounded by a maximum value 
ftmax- Furthermore, the transmission power for each secondary user may be required to be below a 
predefined power limit P max . By incorporating these additional individual constraints imposed on 
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each secondary user, we can model the optimization problem as 

min V e (n,g) 
s.t. C(K,g) < C, k g Z™, g G n%, 

K ^ K max l, ^gf < P max - (17) 

To better understand the optimal resource allocation for cooperative spectrum sensing, we consider the 
following two cases in Scenarios A and B as illustrated in Table ITfl joint optimization of k and g; and 
optimization of either k or g. 



TABLE II 

Optimization Problems for Cooperative Spectrum Sensing 





Case I 


Case II 


Scenario A 


joint optimization of k and g with 
system level cost constraint 


optimization of either k or g with 
system level cost constraint 


Scenario B 


joint optimization of k and g with 
system level and individual constraints 


optimization of either k or g with 
system level and individual constraints 



A. Case I: Joint Optimization of k and g 

1 ) Scenario A: In this case, we consider the optimization in (fT6l) over both n and g. We note that 
(fT6l) is a mixed integer nonlinear optimization problem (MINLP). In general, there is no polynomial- 
time algorithm for solving general MINLPs lfT31 . A potentially clearer insight into the solutions can be 
obtained by considering a convex relaxation for this optimization problem, where we simply relaxed the 
integer constraint of the number of samples: 

min P e (n,g) 

s.t. C(n,g) <C, KeH\, g <Ell\. (18) 

As shown in the Appendix A, (fT8l ) is a convex problem. Thus, it can be solved efficiently using interior- 
point methods or other iterative methods ifToll . This will be a recurring theme in the optimization problems 
we consider in the sequel. In the numerical results, we shall see that the approximation as detailed below 
results in near optimal performance without the curse of complexity. Given this convex optimization 
problem, first we introduce the following lemma. 
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Lemma IV. 1. Optimal solution of (n, g) in (jTS]) should satisfy either 1) Ki > and gi > 0, or 2) ki = 
and gi = for secondary user i. 

Proof: Please see the Appendix B. ■ 
This lemma is not surprising because when one secondary user does not collect the energy samples, it 
will not have anything to transmit to the fusion center. Similarly, when one secondary user decides not 
to transmit the data to the fusion center, it is reasonable to expect that this secondary user should remain 
inactive and not collect local energy samples. Using Lemma II V. 1 1 the optimal solution of («,</) can be 
found as stated in the following theorem. 

Theorem IV.2. Consider the optimization problem in f liffl ), let us define pi = jz y an d assume 

Pi > P2 > • • • > Pn- Then, the optimal solution of («, g) is 

IMC 



Jopt.2) 
K p,i 



(opt,2) 



9 P)l 



0, i>l, 

1/2 



i = 1 



0, i > 1. 



Proof: Please see the Appendix C. 
Given the optimal solution of (n,g), we see that the optimal error probability in (fT8T ) is 

Vc f ji\hi\ 



P^' 2 ) = Q I — max 



2 IcTuVli + \hi\y/co 

Since ([TSj is the relaxation of the MINLP (O, we see that p( opt,1) > p(°P 1 ' 2 ) |J5]. In practice, we may 
consider a floor operation for the number of samples as a suboptimal solution for (fT6l ). i.e., 



„(°Pt,2) 



and ^ b) = <£f' 2) , Vi. (20) 



Let us denote the resulting error probability as pi s ' ub ' ) . Then we see that p(° pt ' 2) < P^ 1 ' 1 ) < P^ sub) . 
Furthermore, when Kp\ is large, based on the first-order Taylor series, we have 

p ( sub ) _ p (o P t,2) = p e(K _ AK , ff )-P e ( K , ff ) « ^^lexp(-5 2 /8)(K 1 + 5 1 )- 2 -> + , 



where 5o = 5 , i7i|/ii|/o' 1 , and 5i = 5f 2 |/ii| 2 /5- 2 . With small value of Aki (normally A«i < 1), it is 
interesting to note that our rounding algorithm is near optimal with large system level cost constraint. 
When C is relatively small, as we will show in our simulations, our proposed suboptimal algorithm can 
also provide a good approximation to the optimal solution. 
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Based on (1201 ). when we jointly design the number of samples and amplifier gains subject to the system 
level cost constraint, only one secondary user needs to be active in the cognitive radio network, i.e., 
collecting local energy samples and transmitting the energy statistic to the fusion center. It is interesting 
to note that this strategy is similar to multiuser diversity where the base station selects the user with the 
highest channel to achieve maximum sum rate capacity ifTTl . In this case, the fusion center will select 
the secondary user with the largest pj to perform local spectrum sensing and data forwarding. This will 
significantly reduce the bandwidth cost for data forwarding. 

Remark: We note that the result in d20l can be implemented in a distributed fashion. The idea is based 
on opportunistic carrier sensing |[T8l or opportunistic relaying |[T9l in which a backoff timer is set to be 
a decreasing function of channel state information. In particular, at the beginning of each sensing time 
slot, the fusion center broadcasts a beacon signal to synchronize all secondary users in the cognitive radio 
network. After estimating the channel gam_| \hi\, the secondary user calculates the control parameter pi 
based on its local received SNR 7, and then maps pi to a backoff timer f(pi) (equal to c/p-i in fT9l , 



where c is a constant). Under a collision free situation, the secondary user with 



argest pi will expire first 
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Note that in this case, 



and perform local energy calculation and data forwarding during this time slot 
fusion center does not need to broadcast the optimal design parameter for each secondary user and this 
will reduce the cooperative sensing cost for broadcasting and reporting. 

2) Scenario B: We examine the optimization (fTTT ) over both k and g. Similar to Scenario A, we first 
consider the relaxation to the original MINLP in (fTTT ). i.e., 

min P e (rc,flO 

s.t. c(K,g) <c, Ken%, gen\, 

K ^ K maxl) £,i9i — T^max- (21) 

Again, we see that this is a convex optimization problem and can be solved by standard methods. Let 
us denote the optimal solution in (fJU) as , g^f^) ■ Similarly, we note that 

Lemma IV.3. Optimal solution of (k,<?) in rt2iD should satisfy either 1) k,i > and gi > 0, or 2) Ki = 
and gi = for secondary user i. 

The proof is similar to that of Lemma II V. 1 1 and thus omitted. With the additional constraints imposed 



3 We assume reciprocity of the uplink and downlink channels between the fusion center and secondary users 1201 . 
Detailed analysis on how to reduce the collision probability for this scheme can be found in (18). 
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on re and g, we see that in general it is difficult to obtain the closed-form solutions for (re, g). Since the 
optimal solution of (re, g) needs to be equal to or greater than simultaneously, we propose a heuristic 
suboptimal algorithm for Scenario B. Specifically, first we assign K max and V max to the secondary user 
with largest pi. If there are remaining resources, we assign n max and V max to the secondary user with 
second largest pi and so on until K max and P max cannot be assigned to any one secondary user. In this 
case, we merely utilize the near-optimal solution in (l20l) to allocate («i,gi) to the secondary user with 
the next largest pi and m = 0, gi = to the rest of the secondary users. Let us denote the suboptimal 
solution as {K^f\ 9pi)- The detailed algorithm for Scenario B is illustrated in Algorithm 1. 

Algorithm 1 Heuristic Suboptimal Algorithm 
Sort pi in a decreasing order. 

for i = 1 to n do 

if c Kmax + ^ max < C then 

C i C C[)K max T'max) 4 ^max> 9i ^ y/'Pmax./ 

else 

Compute Ki and gi from (f20b ; 



Adjust and truncate and gi to guarantee /tj G (0, K max ] and G ^0, \/V max /S,i 
end if 
end for 



and stop. 



B. Cfl^e //: Optimization of Either g or re 

In some applications, either g or re may be fixed for secondary users. For example, local energy 
calculation may be scheduled in a fixed time slot and each secondary user is assigned same number of 
samples. In this case, we need to optimize the amplifier gain to achieve the desired error probability. On 
the other hand, we may need to choose appropriate number of samples when the amplifier gains are fixed. 
Here, we first assume fixed number of samples, i.e., re = re, then we need to minimize the error probability 
by choosing appropriate g. Let us define global transmission power constraint as V to t = C — col T re. We 
now examine both these cases. 

1 ) Scenario A: Here, we minimize global error probability assuming the global transmit power con- 
straint is given as Vtot- We define Zi = gf, ai = k^f and bi = kia1/\hi\ 2 . Then, the optimization problem 
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in (fT6l ) is equivalent to 



mm 



E 



dibs 



i=l 



Zi + 6j 

s.t. £ T 2 < P tot , z y 0. (22) 

It is easy to see that (l22l is a convex optimization problem. After some manipulations, we see that 
the Karush-Kuhn-Tucker (KKT) conditions can be given as 

T 0^ + Ui -\^ i = (23) 

Ao^-T'tot) = (24) 

mzi = 0. (25) 



where Ao > and u; t > are Lagrangian multipliers. First we assume that Ao > and ui = 0, then from 
(l23l) . we see that 



\/a,ibi/{£i\o) - h 



(26) 



where [x] + = max{0, x}. Plugging this into (l24l . we have \/Ao = w 2 i7r——rr- where Sq = {ilzi > 0}. 



Then, we need to determine the set Sq to obtain the closed-form solution for z. To do this, let us define 



Pi = y/bid/ a i- Without loss of generality, we assume /?i < @2 < • • • < n - After some derivations, as 
outlined in Appendix D, we have 

■•• ,i S \f(is)<l,f (is + 1)>1}, f(n)>l 
• • • , n}, otherwise, 



(27) 



where 



Ptot+ELi 6 ^ ' 



Thus, plugging Aq into (1261) . the optimal amplifier gains can be obtained as 



(opt) 

Tp,i 



0, 



1/2 



i e s 

i ^ <S , 



(28) 



(29) 



where j? - 

' E ie s KiVf^Ti/l^il 

Remark: The optimal amplifier gains follow the water-filling strategy, i.e., with larger the chance 
for the secondary user to be inactive is higher, where is a measure of the observation and fusion 
channel quality. Note that oc l/(ji\hi\). Hence, when the local received SNR is low or the fusion 
channel quality is poor, the secondary user tends not to transmit the local calculated energy to the fusion 
center. 



15 



For comparison, we consider two suboptimal solutions for this optimization problem: 1) A simple 
solution is to choose equal transmission power for each secondary user, i.e., g^ = \/V to t/ (n£i); 2) 
Using the Cauchy-Schwarz inequality, we see that P e (Koo) in ® can be minimized when g$ = c^f\hi\ 2 /^, 
where c is a constant. Based on this, we propose an alternate suboptimal solution for amplifier gains, 

(sub) ( ^ I h I ^ \ "^"/^ 

i.e., = I ^„ 7i 7 ' 2 | fe |2^. Vtot ) ■ Let us denote the asymptotic detection probability when hi — > oo 

for these three solutions of amplifier gains as Pi opt ^(Koo), ^^{k-oo) and P^^Koo). Then, we note that 

Lemma IV.4. When f3 2 > fa, P^ equ) ( Koo ) > Pi sub) ( Koo ) > Pi° pt) ( Koo ). 

Proof: Please see the Appendix E. ■ 
2) Scenario B: Next, we minimize global error probability assuming the global transmit power con- 
straint V tot and the individual transmit power limit V ma , x . In this scenario, the optimization problem in 
(fT71) becomes 

a A 



mm 



, Zi+bi 
i=i 



s.t. < Plot, zhO, < P max . (30) 

With the additional constraint in (l30l as compared to (1221 ). the updated KKT conditions are 



+ Ui - Vi£i - A & = (31) 



Vi(&Zi - P max ) = 0, (32) 

where V{ > are Lagrangian multipliers. First we assume that Aq > and Ui — V{ — 0, then from ( 13 1 1 ), 



we see that Z{ = \/ai&i/(£iAo) — bi. Thus, based on the value of \/Ao> we can determine the optimal 
solution of Zi as 

[ 0, if v% > y/di/ibiti) 

^max/Ci, if < VA7 < VaAZi/CPmax + h£i) 



Zi 



\J dih/ (& Ao) - h, otherwise. 
Let us define two disjoint sets for secondary users as Si = {i\zi = 'Pmax/^i} and £2 = < Zi < 
V ma , x /Ci}- Plugging Zi into (f24]), we have 

l^il^max + (1/aAo) Eie5 2 V^iMi ~ T,i<=S 2 b & = ^tot, 



which implies that v/An = ~ — , V";~ 2 _ 

V V U Plot- |5l |P mD +L,(:5 2 b i?i 

In order to determine Si, S2 and \/Ao and thus obtain the closed-form solution for Zi, we propose a 
two-stage generalized water-filling algorithm as follows: 
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1) In the first stage, we aim to determine the set <Si. To do this, let us define = , Without 

loss of generality, we assume fix < P2 < ■ ■ ■ < fin- Then, similar to Scenario A, S\ can be obtained 
by (EH) with 

Ptot 



(33) 



where 5j = {m|/3, m < Pi, i < m < n}. For an outline, please see Appendix F. After S\ is 
determined, we have Zi = V mSLX /£,i, Vi G <Si. 
2) In the second stage, we follow the similar procedure in Scenario A to obtain £2 and Zi for i S±. 
The solution is given in d29l ), except that V tot and n are replaced by V tot — \Si\V mSLX and n — \S±\, 
respectively. 

To summarize, the detailed generalized water-filling algorithm for Scenario B is illustrated in Algorithm 
2. With amplifier gains fixed, we need to optimize the number of samples to achieve the desired error 
probability. In this case, the solutions of the number of samples are similar to those of the amplifier gains 
in both scenarios (with additional relaxation consideration), thus omitted from this paper. 

V. Optimization: Minimization of System Level Cost 

In the section, we aim to minimize the system level cost of cooperative spectrum sensing to achieve a 
targeted error probability. Similar to the optimization problem in Section JV] we consider two scenarios 
which depend on whether additional constraints are imposed or not. For instance, in Scenario A, the 
optimization problem can be formulated as: 

min C{n,g) 

s.t. P e (K,g) < P e , kg ZI, g G U n + , (34) 

where P e is a predefined error probability threshold. Similar to the analysis in Section ITV-All we consider 
the relaxation, i.e., n G 7£™ to this MINLP, and the optimal solution of this relaxation problem is stated 
as follows: 

Theorem V.l. Consider the optimization problem in f li4D and pi as defined in Theorem \IV.2\ Then, 

.(opt) _ 



K d,i 



(opt) 
9d,i 



^v 1 + VSfe)< 



0, i > 1, 



( 1 1 /ITiM 

7?IM 2 V V & ^ 



1/2 

, i = l 

(35) 

0, i>l, 



17 



Algorithm 2 Generalized Water-filing Algorithm 
Stage 1: Sort /3j in an increasing order. 

for i = 1 to do 

Compute f(i) from 031 ); 

if f(i) > 1 then 

Set Si = {1, • • • , i} and stop. 

end if 
end for 
for i G <Si do 

end for 

Stage 2: For j £ Si, sort /3j in an increasing order and set V to t Vtat ~ \Si\V m ax an d n ^— n — |<Si 
for j = 1 to n do 

Compute f(j) from d28l) ; 

if /(j) > 1 then 

Set ^2 = {1, • • • , j} and stop. 

end if 
end for 
for j £ S 2 do 

Compute 7] and Zj from ( |29l ). 
end for 



where e = 4[Q" 1 (P e )] 2 . 

The proof is similar to that of Theorem IIV.2I and thus omitted. Similarly, we may consider a ceiling 
operation for the number of samples as a near-optimal solution for (l34b . Additionally, we see that only 
one secondary user needs be active for collecting the samples for local energy calculation and transmitting 
energy statistics to fusion center. We have separately examined the optimization problem for the remaining 
cases considered in Section[Tvl i.e., when jointly designing k and g for Scenario B; and when designing 
either k or g for both Scenarios A and Scenario B. Due to space limitations, we omit the discussions in 
the paper. 
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VI. Simulation Results 

In this section, we present numerical results for system level performance evaluation and optimal 
design for cooperative spectrum sensing in cognitive radio networks. In the following results, we assume 
°n = a v = 1 an d c o = 1- 

A. Average Error Probability 

In Fig. |3l we plot the average error probability versus the (equal) amplifier gain for all three channel 
scenarios from Table U We see that in the low and moderate fusion SNR regimes, Channel Environment 
II (Rayleigh fading observation channels and AWGN fusion channels) provides the lowest average error 
probability among all three scenarios. Thus, to maintain a desired detection performance, the fusion 
channels need to be as reliable as possible, while the local received SNRs can be dynamic and be used 
to exploit spatial diversity. 

10° 



g. 10-1 

Z5 
co 
.Q 
o 

Q. 
O 

CD 



°- 10- 3 



io- 
ci 5 10 15 20 25 30 

g: amplifer gain (dB) 

Fig. 3. Average error probability for cooperative spectrum sensing. In the simulation, we choose 7 = — 8dB, k — 100 and 
n = 15. 




B. Minimization of Error Probability 

For the optimal system design, we assume n = 6, h = [1.56, 1.99, 0.37, 1.52, 0.39, 1.98] T and 7 = 
[-8.86, -15.23, -7.21, -5.09, -10.00, -10.97] 1 (dB). Here, we define the global fusion SNR as SNR = 
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Vtot/o'v- For comparison, we consider equal number of samples and amplifier gains as a suboptimal 
solution. 

In Fig.Hl we plot the error probability versus system level cost constraint in Case I for joint optimization 
of k and g for both Scenario A and B. In this simulation, we utilize standard MINLP methods 11211 for 
optimization problem in (fT6l ); the closed-form solution {K^f' 2 \ 9^ ^) m Theorem IIV.2I for the convex 
relaxation of the optimization problem in (fT8T ); and our proposed suboptimal solution («p™ b \ Sp™ ) in (l2Ql > 
in Scenario A and interior-point method to solve the optimization problem in Scenario B. As expected, 
we see that in Scenario A, the error probability of optimization problem in ( fT6l ) and its relaxation in (fT8T ) 
converges, even with relatively small system level cost constraint. Also, our proposed suboptimal solution 
in (l20l ) is near optimal as previously mentioned. Furthermore, we observe that the error performance is 
degraded with the additional constraints in Scenario B. Additionally, our proposed suboptimal algorithm 
in Scenario B has negligible performance loss compared to the optimal solution. 




■| Q- 4 I I I I I I I I I I 

14 16 18 20 22 24 26 28 30 32 34 

C ■ system level cost constraint (dB) 



Fig. 4. Case I: error probability for different solutions of (n, g). In Scenario B, we choose K max = 0.2[C/coJ and P max = 0.2C. 

Fig.[5]shows the error probability versus total number of samples in Case II (optimization of g given k). 
As expected, we see that the optimal solution provides superior performance to suboptimal solutions. From 
the plots, we also observe that with additional individual constraints, the optimal solution for Scenario 
B performs worse than that of Scenario A. Furthermore, when total number of samples increases, we 
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see that the error probability approaches the asymptotic bound. In particular, Pe («oo) > Pe Sub ^(Koo) > 
Pi° pt ^(Koo) as stated in Lemma Irv74] 




Fig. 5. Case II: error probability for different solutions of g. In the simulation, we choose SNR = 25dB and fixed number of 
samples Hi = \_Ktot/n]. In Scenario B, we choose "Pmax = 0.4"P to t. 



VII. Conclusions 

In this paper, we present the performance evaluation and optimal design for spectrum sensing in the 
cognitive radio networks. We first analyze the average error probability by considering a range of channel 
realizations between the primary user and the secondary users and between the secondary users and the 
fusion center. Then, we investigate the optimization problems for spectrum sensing. In particular, when 
jointly designing the number of samples and amplifier gains, we demonstrate that only one secondary 
user needs be active, i.e., collecting local energy samples and transmitting energy statistic to fusion 
center. Furthermore, we derive closed-form expressions for optimal solutions and propose a generalized 
water-filling algorithm when number of samples or amplifier gains are fixed and additional constraints 
are imposed. 
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VIII. Appendix 
A. Proof of Convexity of Optimization Problem (1181) 

Proof: Let us define Zj = gf, p { = 5-%/(jf\hi\ 2 ), q { = I/7? and T^ki, zi) = p .£+ q . z . ■ To simplify 

our analysis, when m = z% = 0, we assume ^(/tj, 2j) = Q3 Then, the optimization problem ( fT8l ) becomes 

n 

max V" Ti{Ki,Zi) 

K,Z * ' 

8=1 

s.t. c 1 t k + £ t ,z < C, K ^ 0, 2 ^ 0. (36) 
After some manipulations, we see that the Hessian of Fi{Ki, z\) is given as 

2pi% 



V Ti{Ki,Zi) 













Zi 


Ki 




Ki 



-< 0. 



(PiKi + qiZi) 3 

Thus, Ti(Ki, Zi) is a concave function, which indicates that the objective function in (l36l ) is also concave. 
This completes the proof. ■ 

B. Proof of Lemma \1V.1\ 

Proof: We prove this lemma by contradiction. First we assume that (k, z) with Ki = 0, Zi > or 
Ki > 0, Zi = for secondary user i is the optimal solution for (l36l ). Let us define the optimal value is 
p*. Since KjZj = 0, the objective function remains unchanged in (l36l ). Then, the optimization problem 
becomes 

max :,;,,-^i'0-"-j) 

Kj > 0, > 0, Vj 7^ i. (37) 

where C = C — when k, = 0, Zi > 0, or C' = C — co^i when k» > 0, = 0. In either case, we 
see that C' < C. To prove this lemma, we need to find a substitute solution (k',z') with optimal value 
p'* > p*. To do this, let us replace the solution for secondary user i as k\ = z[ = 0. In this case, the 
optimization problem becomes 

s.t. c E"=ij#i K i + £"=ij#i < c 

> 0, > 0, Vj / i. (38) 



5 In practice, this assumption can be alleviated by adding a sufficiently small constant in the denominator. 
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Then, we see that it is equivalent to proving that the optimal value p'* in (l38l ) is greater than p* in (I37T ). 
Since the objective and constraint functions in these two optimization problems are identical, this can be 
easily proved by convex relaxation in optimization problem, which implies that we can find a substitute 
solution (k',z'), i.e., p'* > p*. This contradicts the assumption that (k,z) is the optimal solution and 
we can conclude the proof. ■ 

C. Proof of Theorem \IV.2\ 

Proof: The Lagrangian function of (l36l ) can be given as 

n 

C(K, Z, Xq,U,v) = — ^2 \~ \q(cq\ T K + £ T z) — U T K — V T Z — XqC, 

where Ao > 0, Uj > and V{ > are Lagrangian multipliers. Here the KKT conditions are 

(^fep + «* ~ ^Ao = (40) 
X (c I t k + $, t z-C) = (41) 
UiKi = 0, ViZi = 0. (42) 

From Lemma IIV.1[ we see that U{ and v-i need to be or greater than simultaneously. First we assume 
Ui = Vi = and Ao > 0, which indicates that m > and Zi > 0. Then from ( f39b and (1451 . we have 
= UiKi, where = \/coPi/{qi£i). Plugging this into (f36l >. the original optimization problem becomes 

max Yjiex s u K i 

S - L Y,i&T S 2i K i ^ Kj > 0, Vi G X, (43) 

where X = > 0,Zj > 0}, sij = and S2i = Co + ^cjj. Since adding zero will not 

change the objective function and constraints in (l43l) . we can rewrite (l43l) as 



max s^k 



s.t. s^k < C, k ^ 0. (44) 

This is a classic linear optimization problem; thus we can solve this easily. Since the vertices of the 
polyhedron are the basic feasible solution for linear optimization problem E2l . the optimal solution of 
(l44l) suggests that only one of m is non-zero while others are all zero. Let us define pi = su/s2i and 
assume p\ > p2 > • • • > p n . Then, the optimal solution of can be given in Theorem IIV.2I This 

completes the proof. ■ 
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D. Solution for Set Sq 

Here we follow the analysis in |23l to find Sq. From d26l ), we see that in order to guarantee %i > 0, 



we need to have ^/Ao < V a i I » which indicates /(i) < 1 for some is. Then, the problem can be 
stated as: given fa < fa < ■ ■ ■ < /3 n , /fe) < 1 and /(i^ + 1) > 1, we have 

1) f(i) is an increasing function of i for i < is; 

2) f(i) > 1 for i > i 5 . 

Proof: It is straightforward to show that /(l) < 1. This indicates that So ^ and thus there exist 
feasible solutions for z. When i > 1, we have 

71 + ; ~~ £}=i &i&+Pw+&t+i&+i 
> ft Ej = i \/ajbjgj+6i+igi+i 

(») | /(i), i < is 
[ 1, i > «5- 

The first inequality in (a) is valid since when x/y < 1, we have (x + c) / (y + c) > x/y, where x,y,c> 0. 
Then, we see that f(i) is an increasing function of i for i < is- The second inequality in (a) is valid 
since when x/y > 1, we have (x + c)/{y + c) > 1. This indicates that when f(i) > 1, f(i + 1) > 1 for 
i > is- This completes the proof. ■ 

E. Proof of Lemma \1V.4\ 

Proof: From Section HV-B II we see that 

f(2) = bib. + h& + (fa ~ fhWathtj 

As — > oo, we have ai,&i — > oo. This implies (/% — > With fa > fa, /(2) > 1 

and S = {1}. Then, Pi opt) (/c°o) = Q (^{V to t maxjfl,}) 1 / 2 ), where ^ = 7 ?N7&- Furthermore, 
Pi Sub) (^oo) = Q (^(PtotllflllVa^)) 172 ) andPi equ) ( Koo ) = Q (^(T^l^/n) 1 ^). Since max{^}- 
(1 T #) > ||#|| 2 and n||0|| 2 > (1 T 0) 2 , we can conclude the proof. ■ 

F. Solution for Set S\ 

Similar to the solution for Sq, we need to show that: given fa < fa < • • • < /fo) < 1 an d 
/fe + we have 

Property F.l: f(i) < 1 for i < 2,5; 
Property F.l: f(i) > 1 for i 5 < i < [^-J . 
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Proof: To prove Property F.l, we consider 4 cases which depend on the values of ft and ft: 1) 
Si = <Si_i U 5,' \ {i}, 2) Si = U S[, 3) 5i = Si-x \ {i}, 4) Si = Si-i, where = {m|ft_i < ftn < 
ft, i < m < n}. Now we start with case 1). In case 1), we have ft < ft_i and S' { ^ 0. Furthermore, 
we note that 

A— 1 XymeSi-i V a mbm£,m 

"msri 



< ft£ 

The last inequality is valid since when m E 5,', /3 m < ft, we have 

mSiS,' mdSl m ScS l ' 

After some manipulations, when /(i) < 1, Vz < i^, 



?/. n s ^ ft E m6lSl Vam&m^m + C1 

A* -1 ) - 7r3 ^5 a c \ ■ < lj 

where ci = (V max + — J2 m ^s' bm£,m- The last inequality is valid because when x/y < 1, we have 
(x + ci)/(y + ci) < 1, where x,y > and ci > — x. Similarly, we see that for other three cases, we 
also have f(i — 1) < 1. 

Now let us prove Property F.2. Similar to Property Fl, we have 

r , . ft Eme^ VambmZm ~ C 2 

f( l + 1 ) ^ 7r3 ^5 — r~r~\ - L 

(>tot - l^max + Z^meS, b mt,m) ~ C 2 

where c 2 = (P ma x + ^i+iCi+i) ~~ SmeS' b m ^ m . The last inequality is valid because when x/y > 1, we 
have (x — c 2 )/(y — c 2 ) > 1, where x,y > and c 2 < y. Similarly, we see that for other three cases, we 
have f(i + 1) > 1. This completes the proof. ■ 

References 

[1] S. Haykin. Cognitive Radio: Brain-Empowered Wireless Communications. IEEE J. Sel. Areas Commun., 23(2):201-220, 
Feb. 2005. 

[2] M. Nekovee. Current Trends in Regulation of Secondary Access to TV White Spaces Using Cognitive Radio. GLOBECOM 
2011, pages 1-6, Dec. 2011. 

[3] Z. Quan, S. Cui, and A. H. Sayed. Optimal Linear Cooperation for Spectrum Sensing in Cognitive Radio Network. IEEE 

Journal of Selected Topics in Signal Processing, 2(1):28— 40, Feb. 2008. 
[4] J. Ma, Y. G. Li, and B. H. Juang. Signal Processing in Cognitive Radio. Proc. IEEE, 97(5):805-823, May 2009. 



25 



[5] A. Ghasemi and E. Sousa. Collaborative Spectrum Sensing for Opportunistic Access in Fading Environments. DySPAN 

2005, pages 131-136, Nov. 2005. 
[6] G. Ganesan and Y. G. Li. Cooperative Spectrum Sensing in Cognitive Radio - Part I: Two User Networks. IEEE Trans. 

Wireless Commun., 6:2204-2213, June 2007. 
[7] J. Unnikrishnan and V. V. Veeravalli. Cooperative Sensing for Primary Detection in Cognitive Radio. IEEE Journal of 

Selected Topics in Signal Processing, 2(1): 18—27, Feb. 2008. 
[8] S. Zhang, T. Wu, and V. K. N. Lau. A Low-Overhead Energy Detection Based Cooperative Sensing Protocol for Cognitive 

Radio Systems. IEEE Trans. Wireless Commun., 8(1 1):5575-5561, Nov. 2009. 
[9] S. Appadwedula, V. V. Veeravalli, and D. L. Jones. Energy-Efficient Detection in Sensor Networks. IEEE J. Sel. Areas 
Commun., 23(4):693-702, April 2004. 
[10] Z. Quan, S. Cui, A. H. Sayed, and H. V. Poor. Optimal Multiband Joint Detection for Spectrum Sensing in Cognitive 

Radio Networks. IEEE Trans. Signal Process., 57(3): 1128-1 140, Mar. 2009. 
[11] A. Sahai, N. Hoven, and R. Tandra. Some Fundamental Limits on Cognitive Radio. Proc. Allerton Conf. Communication, 

Control, and Computing, pages 131-136, Oct. 2004. 
[12] A. T. Hoang, Y. C. Liang, D. Wong, Y. Zeng, and R. Zhang. Opportunistic Spectrum Access for Energy-Constrained 

Cognitive Radios. IEEE Trans. Wireless Commun., 8(3): 1206-121 1 , March 2009. 
[13] M. Alouini and A. J. Goldsmith. A Unified Approach for Calculating Error Rates of Linearly Modulated Signals over 

Generalized Fading Channels. IEEE Trans. Commun., 47(9): 1324-1334, Sept. 1999. 
[14] I. S. Gradshteyn and I. M. Ryzhik. Table of Integrals, Series, and Products, 6th ed. San Diego, CA: Academic Press, 
2007. 

[15] L. A. Wolsey. Integer Programming. Wiley, 1998. 

[16] S. Boyd and L. Vandenberghe. Convex Optimization. Cambridge, U.K.: Cambridge Univ. Press, 2003. 
[17] R. Knopp and P. Humblet. Information Capacity and Power Control in Single Cell Multiuser Communications. ICC, pages 
331-335, June 1995. 

[18] Q. Zhao and L. Tong. Opportunistic Carrier Sensing for Energy Efficient Information Retrieval in Sensor Networks. 

EURASIP J. Wireless Commun. Netw., 2:231-241, April. 2005. 
[19] A. Bletsas, A. Khisti, D. P. Reed, and A. Lippman. A Simple Cooperative Diversity Method Based on Network Path 

Selection. IEEE J. Sel. Areas Commun., 24(3):659-672, March 2006. 
[20] T. S. Rappaport. Wireless Communications: Principles and Practice. Prentice Hall, Upper Saddle River, NJ, USA, 1996. 
[21] I. E. Grossmann and Z. Rravanja. Mixed-Integer Nonlinear Programming: A Survey of Algorithms and Applications . 1997. 
[22] D. Bertsimas and J. Tsitsiklis. Introduction to Linear Optimization. Athena Scientific, Belmont, Massachusetts, 1997. 
[23] S. Cui, J. Xiao, A. J. Goldsmith, Z.-Q. Luo, and H. V. Poor. Estimation Diversity and Energy Efficiency in Distributed 

Sensing. IEEE Trans. Signal Process., 55(9):4683-4695, Sept. 2007. 



